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1. Introduction 



To be able to predict the evolution in time for a Hamiltonian system is of central importance 
to most practical physical problems. For the standard situation, i.e. when the Hamilto- 
nian is Hermitian, there exist well developed frameworks. Having for instance in mind to 
anticipate the response of an atomic system described by a Hermitian Hamiltonian when 
5^ \ it is subjected to an external time-dependent laser field is an intensively studied problem 

in the weak and recently also in the strong field regime. In the former case Fermi's golden 
rule is for instance one of the central results, whereas the latter case leads to interesting 
new phenomena such as high harmonic generation [Q, above threshold ionization [^ and 
stabilization |^ ^, |5| , see also |6|, for strong field phenomena in general. 

Less developed is the situation regarding non-Hermitian Hamiltonians. Depending on 
the nature of their eigenvalues non-Hermitian Hamiltonian systems can be investigated 
in various fundamentally different ways. When the corresponding energy eigenvalues are 
complex one may essentially keep the standard framework and accept the fact that the non- 
Hermitian nature of the Hamiltonian will lead to decaying states and wavef unctions. Vari- 
ous investigations concentrate on that particular setting |@, |^, IC, 11, ^, 13, 14|. However, 



more recently it was observed that a large class of non-Hermitian Hamiltonians possess real 



and positive spectra |15|. For that situation it is natural to demand preservation of prob- 



ability density, which is not guaranteed even when the Hamiltonian is time- independent. 
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By now there exists a considerable amount of results regarding this situation for various 
single particle systems |jl^, ^, |l^, |l^, 20 1 quantum field theories^ [23, 24, |^] and also 



integrable many particle systems ||2^, ^] . So far most effort has gone into the development 
of a proper quantum mechanical framework for such systems. The main purpose of this 
paper is to extend these treatments, such that they include a proper description for the 
time-evolution for non-Hermitian Hamiltonian systems. 

Our manuscript is organized as follows: In order to establish our notation and to 
highlight the key concepts we review in section 2 the main characteristics for a consistent 
quantum mechanical description involving non-Hermitian Hamiltonians. In section 3 we 
generalize the scheme to construct time-independent pseudo Hermitian Hamiltonian sys- 
tems and provide a systematic procedure, which leads to closed formulae involving Euler's 
numbers for the equivalence pairs of Hermitian and non-Hermitian Hamiltonians, h and 
H, respectively. Subsequently we employ that scheme to compute various new equiva- 
lence pairs in its exact and perturbative form needed afterwards for the time-dependent 
treatment. In particular, we generalize non-Hermitian perturbations of the harmonic oscil- 
lator to anharmonic oscillators with a wider class of perturbations in an exact formulation. 
Amongst those new non-Hermitian Hamiltonians is a doubly graded generalization of the 
Swanson Hamiltonian. Perturbatively we also compute Hermitian counterparts for the 
harmonic oscillator with an additional igx'^-ieiTa for generic n, hitherto only studied for 
specific cases. In particular, for the case n = 3 we provide the explicit formula for all 
wavefunctions up to order g^^ which turn out to be far simpler than their non-Hermitian 
counterparts. In section 4 we discuss the time evolution for non-Hermitian Hamiltonian 
systems in various different gauges and investigate a time-dependent perturbation theory. 
As a particular example we employ the formalism to compute some transition probabili- 
ties for the harmonic oscillator perturbed by igx^ in an external laser field. We state our 
conclusions in section 5. 

2. Quantum Mechanics involving non-Hermitian Hamiltonians 

The possibility that non-Hermitian Hamiltonian systems can possess discrete eigenstates 
with real positive energies has already been indicated by von Neumann and Wigner |^8| 
almost eighty years ago. More recently this type of systems are under more intense scrutiny 
and nowadays the properties of these so-called BICs (bound states in the continuum) are 
fairly well understood for many concrete examples p^, pH, pl| together with their bi- 
orthonormal eigenstates p2| , p3[ . 

Whereas the above type of Hamiltonians only possess single states with these "strange 



properties" [^8[, it was observed eight years ago by Bender and Boettcher [15| that Hamil- 
tonians with potential terms V = x^{ixY for v > Q possess an entirely real and positive 
spectrum. Since that discovery non-Hermitian Hamiltonians, in the sense that W ^ H, 
are under intense investigation. Initially the reality of the spectrum was attributed to 
the PT-symmetry of the Hamiltonian. In fact, when the wavefunctions are simultaneous 



^In 1+1 dimensional quantum field theories non-Hermitian Hamiltonian systems are known to be mean- 
ingful for some time pil, [22|. 
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eigenstates of the Hamiltonian and the PT-operator one can easily argue that the spec- 



trum has to be real [34|. However, despite the fact that [PT,H] = 0, this is not always 
guaranteed as the PT-operator is an anti-linear operator p^ . As a consequence one may 
also encounter conjugate pairs of eigenvalues for broken PT-symmetry [^. To determine 
whether the PT-symmetry is broken or not one may use various techniques to verify this 
case- by-case ||T^, 35 1. The central problem arising in this context is that inner products 



of wavefunctions constituting solutions of the time-independent Schrodinger equation in- 
volving non-Hermitian Hamiltonians become indefinite, which is due to the fact that the 



wavefunctions have to be simultaneous eigenfunctions of H and the PT-operator |36|. Ben- 
der, Brody and Jones [g4[ solved this problem consistently by introducing a new type of 
inner product 

m<^')^p^:={CPT\^)f ■\<^'), (2.1) 

which then indeed leads to a positive definite metric, that is ($„| ^rn)cPT ~ ^nm when 
labelling the energies by increasing values e„- This inner product inherits one complication, 
which is already present when solving the eigenvalue problem, namely that eventually the 
wavefunctions $ no longer vanish for |x| — > oo. In that situation one has to integrate 
within wedges bounded by the Stokes lines in the complex x-plane [|l5| . A further initial 
drawback of this formulation was that the C-operator C{x, y) = ^^ <l>„(x)<&„(y) needed to 
be determined dynamically, which requires in principle the knowledge of all wavefunctions. 
Meanwhile also alternative methods have been developed to compute C and this is no 
longer a real obstacle. For instance, noting that C is a symmetry of the Hamiltonian 
and in addition an involution, one may compute it alternatively by solving the algebraic 



equations |37] 



[C, H] = 0, [C, PT] =0 and C^ = 1 . (2.2) 

Even before the discovery of [^] and the introduction of the CPT-inner product ( ^.l| ) 
there have been very general considerations addressing the question of how a consistent 
quantum mechanical framework can be constructed for non-Hermitian Hamiltonian systems 



1 38]. It was understood that quasi-Hermitian (pseudo-Hermitian) systems would lead to 



positive inner products. Subsequently this was further developed by Mostafazadeh |39 



40, 41, ^3], who proposed that instead of considering PT-invariant Hamiltonians one may 



investigate pseudo-Hermitian Hamiltonians satisfying 

h = r]Hr]-^ = h) = r]-~^H^r] ^ H^ = rfHrj-'^, (2.3) 

with if = Tj. Since the Hermitian Hamiltonian h and the non-Hermitian Hamiltonian H 
are related by a similarity transformation, they belong to the same similarity class and 
therefore have the same eigenvalues. The corresponding time-independent Schrodinger 
equations are then simply 

h(l) = ecj) and H^ = e$, (2.4) 

where the wavefunctions are related as 

$ = 7?"V- (2.5) 
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Having real eigenvalues for the Hermitian Hamiltonian h then guarantees by construction a 
positive spectrum also for H. In fact, this formulation is more general than demanding the 
Hamiltonian to be PT-symmetric, which is only a sufficient, but not a necessary condition 
for the spectrum to be real for unbroken PT-symmetry of the wavefunctions. In addition 
the formulation which involves pseudo-Hermitian Hamiltonians is more intuitive as the 
reality of the spectrum of H is completely evident. Inner products for the wavefunctions 
$ related to the non-Hermitian Hamiltonian H may now simply be taken to be 



($!$') :=(<I>|r?2$'>, (2.6) 



where the inner product on the right hand side of ( |2.6D is the conventional inner product 
associated to the Hermitian Hamiltonian h. In case the similarity transformation ( |2.3| ) 
holds, the Hamiltionian H is PT-symmetric and when in addition the solution to ( p.2| ) is 
taken to be C = rj~^P, the CPT-inner product (^.l|), the //-inner product ( ^.6| ) and the 
conventional inner product related to the Hermitian Hamiltonian coincide 

{<^\<^')^p^ = {^\^')^ = {cp\cl>'). (2.7) 

With regard to ( ^.71) one may wonder why one requires the CPT -inner products when 
one may in fact use the T/-inner products, or even more radically why one needs the non- 
Hermitian formulations at all when they can always be related to the standard inner prod- 



ucts. In fact, these issues are quite controversially discussed at present ||4^, ^, ^, ^46| . 
With regard to CPT versus pseudo-Hermiticity, our point of view is that despite the lim- 
ited restrictive power of PT-symmetry, in particular the fact that is does not guarantee 
a positive spectrum, it is a very good guiding principle to select potentially interesting 
non-Hermitian Hamiltonians on the classical level, e.g. for many-particle systems |2(:, 27|. 



This property can be read off directly from a classical Hamiltonian, whereas even when 
one has identified such Hamiltonians, a proper analysis requires the construction of the 
similarity transformation 7] of the CPT-operator, which is usually not evident a priori. 
With regard to the inner products, it appears far easier to construct ry rather than the 
CPT-operator. One apparent virtue of the non-Hermitian formulation, using CPT or 77- 
inner products, is that in this way we may relate simple non-Hermitian Hamiltonians to 
fairly complicated Hermitian Hamiltonians. It is sometimes argued that the computations 



in the non-Hermitian framework are simpler to perform |45|, but this statement has been 
challenged [^] . Certainly, as we will see below, this feature can not be elevated to a general 
principle. We will see that even when the non-Hermitian Hamiltonian looks simpler than 
its Hermitian counterpart, this is not true for the corresponding wavefunctions, which still 
take on a simpler form in the Hermitian formulation. Furthermore, we find here in addition 
that the time-dependent non-Hermitian formulation will always be more complicated or 
at most of equal degree of complexity than the Hermitian one. It appears to us that the 
best strategy is to make use of both worlds and switch to one or the other formulation 
depending on the specific problem at hand. 

The main purpose of this paper is to investigate how these frameworks may be trans- 
lated to the situation when the Hamiltonians become genuinely time-dependent. 
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3. Construction of pseudo Hermitian Hamiltonians 

Accepting that a non-Hermitian formulation of quantum mechanics is more straightforward 
in a pseudo-Hermitian formulation rather than a CPT-scheme, the question with regard 
to ( |2.3| ) then arises of how to construct Hamiltonians h and H belonging to the same 
equivalence class. Supposing that the similarity transformation can be realized by using an 
operator of the form 77 = exp(g/2), the relation (|2.3D implies by standard Baker-Campbell- 
Hausdorff commutation relations that 

1 1 °° 1 

H^ = H+[q,H] + -[qM,H]] + -[qMAa.H]\] + --- = Y.^:'^^^^H). (3.1) 

n=0 

For convenience we have introduced here a more compact notation for the n-fold commu- 
tator of the operator q with some operator x as 



.(") 



(x):=[g,[g,[(7,...[(?,x]...]]]. (3.2) 



Taking now the non-Hermitian Hamiltonian to be of the form H = hQ + ihi, with /iq = h^, 
hi = h\ the relation ( ^ ) acquires the form 

i [Q, ho] + ^ [q, [q, ho]] + ^ [q, [q, [q, ho]]] + ... = 2hi + [q, hi] + - [q, [q, hi]] + ... (3.3) 

In solving this equation one may start from different given quantities. For instance one 
may solve for hi with given ho, q, see section 3.1, or one may solve for q with given ho, hi, 
see section 3.2. We will not treat the remaining possibility to construct /iq for given hi,q. 
Furthermore, we will also not discuss the interesting possibility to exploit the isomor- 
phism between commutator relations and Moyal products [^, 48 1. This relation allows to 



translate the commutator relation into a differential equation for rj, which may be solved 
subsequently. 

3.1 Exact similarity relations 

Some simple exact solutions to ( |3.3| ) can be found easily in a quite systematic way by search- 
ing for vanishing multi-commutators. For instance, if for a given Hermitian Hamiltonian 
ho we can find some q such that its triple commutator with ho is vanishing 

[q,[q,[q,ho]]]=0, (3.4) 

we can define the non-Hermitian part of H as 

1 
hi = -[q,ho], (3.5) 

such that the relation ( |3.3| ) is solved exactly. According to (p]^) the Hermitian counterpart 
h of the non-Hermitian Hamiltonian H is then computed to 

h = r](ho--[q,ho]j r]'^ = ho - - [q, [q, ho]] ■ (3.6) 
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We can generalize this construction procedure to any vanishing n-fold commutator of q with 
Hq. Assuming for this that the sum in ( ^ ) terminates at some stage, i.e. Cq (/lo) = 0; 
we make the following ansatz for the non-Hermitian part of H 

^^ = ^E54"^(^o), (3.7) 

n=l 



where the constants k„ E M are to be determined such that the relation (3^) is solved 



exactly. The substitution of (|3.7D into (|3.3|) then yields 



£ , l i-1 



E ^/^'(ho) = E 2^4")(/.o) + E E ^.^-r-H^^o). (3.8) 



n! ^ ^ — ' n! " ^ — ' ^ — ' nl m\ 

n=l n=l n=l m=l 



Reading off the coefficients of equal n-fold commutators from this equation produces a 
recursive equation for the constants k^ 



-I -I n-l , 

1 1 v-^ / n 



2 2 ^^ \m 

m=\) 



,„_--- E ( K^ for 1 < n < £. (3.9) 



With kq = 0, we can solve ( |3.9D iteratively and find that all coefficients Kn for even n 
vanish, whereas the remaining ones become 

1 11 17 31 691 , , 



Finally we compute from ( |2.3|) the Hermitian counterpart h of H to 

h = r] 
where the constants An are related to the k„ as 



^^0 - E S4'^^(^o)) .-^ = E ^cW(^o), (3.11) 

V n=l / n=0 



A„ = 1-X:2™Qk... (3.12) 

m=0 ^ ^ 



Using the above solutions for the «;„ ( p. 101 ), we find that only coefficients A„ with n even 
are non- vanishing 

Ao = l, A2 = -1, A4 = 5, Ae = -61, As = 1385, Aio = -50521, .. . (3.13) 

In fact we observe that these constants are very closely related to Euler's numbers En as 
A2n = (— 1)"'-E'„ for n = 1, 2, 3, . . .With this identification we may alternatively solve the 
equations (|3.12| ) for the constants Km, such that they are also expressed in terms of Euler's 
numbers 

, [("+l)/2] . X 

m=l ^ ^ 

Here [x] denotes the integer part of x. 
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Thus given some Hamiltonian /iq, which constitutes the Hermitian part of a non- 
Hermitian Hamiltonian H = hg + ihi, together with an operator q satisfying Cq (ho) = 
for some finite integer i, the above procedure provides a systematic way to compute pairs 
of Hamiltonians 

h = ^o + E^-^^-tHho) and H = ho- Y. jS^/t-'^i^o), i^-l^) 

ra=l ^ '' n=l ^ '' 

with h = h'^ and H ^ H^ , which belong to the same similarity class related by the adjoint 
action of r/ = exp(g/2) according to (p.3|). The closed formulae in ( p. 15 ), together with the 
line of arguments leading to them, appear to be new. 



3.2 Perturbative similarity relations 

Often one has a different type of starting point as in the previous subsection and would like 
to construct h for a completely specified non-Hermitian Hamiltonian H, that is for given 



/iQ and hi. In that case we have to solve the commutator relation (3/?) for q and some £. 
Whenever this is not possible in an obvious manner, one can resort to perturbation theory 
as originally proposed by Bender, Brody and Jones [^ (see also |4^). To develop this 
one makes a further assumption on the form of the similarity transformation r] = exp(q/2), 
namely 



£9'"-^2n-i. (3.16) 



n=l 



One may argue, by demanding PT-invariance [37|, that the powers oi g have to be odd. 
Here we want to guarantee pseudo-Hermiticity and therefore present a slightly different 
argumentation. We assume the following dependences on the coupling constant g 

Vi-9)=v{9r\ h{g)=h{-g) and H^g) = H{-g). (3.17) 



The first equation is obviously satisfied by the ansatz ( 3.16| ), whereas the second and third 



are supported by the examples presented below. Using then ( p. 17 ) the pseudo-Hermiticity 



H\g)=r^{gfH{g)r^{gr^ (3.18) 

simply follows from 

h{g) = r,{g)H{g)rj{g)-' = h{-g) = rj{-g)H{-g)r,{-g)-' = r^{g)-^ H\g)i^{g) . (3.19) 



Returning to the discussion of perturbation theory, we see that with the ansatz ( 3.1(i| ) 



the multi-commutator Cg (/iq) will be at least of the order 0{g^). This means that a 
precision of order 0{g ) corresponds to Cq yio) = 0, such that the above arguments 
apply and from (^j) we obtain 

^1 = ^ES5" E 4r^ (4r^ (...4^^)(m)) +o(/+^). (3.20) 

n=l ' n\+n2+ ■ ■ ■+ni=n 
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Solving these equations order by order yields the set of equations 



[/io,gi] = —hi, 
9 



[ho,qb] 



Qg 



Qg 



cff(c«(/^i))+4?(e«(/^i))-^c(t)(/.i) 



(3.21) 
(3.22) 
(3.23) 



which can be used to determine the unknown quantities qi for 'i < i < H. recursively, 



as already noted in |37|. Having determined the qi to the desired order the Hermitian 
counterpart h to H results from ( 3.11| ) to^ 



[£/2] 



""E 



n=0 



4"(2n)! ^ 



E 



^qi V "^132 



(ni) ( Jn2) 



. . .C 



{•"-1)1 
Qi 



{ho)]]+0{g 



i+i\ 



(3.24) 



"i+"2+---+'^£=2ra 



We present here now various time-independent Hamiltonians which belong to the same 
equivalence class and which we discuss below in its time-dependent variant. 



3.3 Non-Hermitian Hamiltonians and their Hermitian counterparts 

We will consider some non-Hermitian perturbations of the harmonic or anharmonic oscil- 
lators depending on a real coupling constant a S M 



^» = ^P^ + f-^ 



(3.25) 



where x and p are operators obeying the standard canonical commutation relation [x,p] = i. 
Throughout this paper we use atomic units h = e = nie = ca = 1. In the last equation a 
is of course the fine structure constant and not the coupling constant in ( 3.25| ). 



3.3.1 Anharmonic oscillator perturbed by iJ29p^^ 

We take as a starting point the harmonic oscillator /iq = h2{a) and q = fip, where n is 
a real constant which needs to be determined. It is then easily checked that the triple 



.(3) 



commutator c\ (/iq) indeed vanishes. Therefore we evaluate from (3.5) 



hi = -afix = gx, 



(3.26) 



where we introduced a new coupling constant 5 G M to simplify the notation. The non- 
Hermitian Hamiltonian is then of the form 



H{a, g) = -p^ + -x^ + igx. 



(3.27) 



^ After completion of this work we were notified by H.F. Jones that he also noted the occurrence of 



Euler's numbers in front of the multi-commutators of the perturbative expressions (3.21)-( |3.23| ) [pO|. Here 
they are a consequence of our general formulae (3.15). 
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According to ( p. 61) we compute next the Hermitian counterpart of H{a, g) to 

h{a,g) = r]H{a,g)r]-^ = \p' + ^x^ + \^. (3.28) 

2 2 2 a 

with 

7] = e^P. (3.29) 

This equivalent system of Hamiltonians (|3.27] ) and ( p.28|) follow also directly from ( 3.15| ) 



and for a = g = 1, they can already be found for instance in |45]. 

We may now easily generalize this system by taking for the Hermitian part of the 
Hamiltonian H the Hamiltonian /iq = h^{oi) as starting point. We compute 

cS'«(«)) 4 '--''"''^^" !" '^'"^". (3.0, 

I U tor m > n 

and thus we can take here £ = n as a cut-off condition in order to compute hi from (p. 7]). 
For fi = 2g we then find 

-I . [("+l)/2] , X 2m-l / X 

where the constants Km are determined from ( p. 141 ) as sums over Euler's numbers. The 
Hermitian counterpart results from ( 3.11| ) or ( |3.15| ) to 



V2] /„ X 2m 



Clearly, since c\l {h^{a)) = for some finite values of k,n,m, we can generalize this and 

take q = Ylrn=i f^mP"^ to construct further conjugate pairs of Hamiltonians. Notice that 
the dependence of h^'^{a,g), H^^{a,g) and r]{g) on the coupling constant g respects the 
aforementioned identities ( |3.17| ). 



3.3.2 Generalizations of the Swanson Hamiltonian 

Next we start again with /iq = h^{a), but change the operator q in the similarity trans- 
formation to qm = /^m^"*- It is easy to observe that in this case Cn\j.m{h^{a)) = for all 
m, n > 0. Therefore when taking fi^^ = 2g/m relation ( p. 5] ) yields 

^n,m(«' 9) = \v^ + f ^" - if (px— 1 + x^-^p), (3.33) 

and consequently (|3.6|) gives 



^n,m(">5) = /i°(a) - - [qm, [gm,/i°(a)]] = vHlrni<^)v ^ (3.34) 

= ^/ + fx" + ^5V-2. (3.35) 
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Notice that when specializing to n = m = 2, we obtain the harmonic oscillator 

hl^{a,g) = hl{a + g^) (3.36) 

and H22 becomes the Swanson Hamiltonian ||5l|, ^ 

H{a,g) = hl{a)-i-{xp + px), (3.37) 

upon changing the conventions for the coupling constants. We further note that only for m 
odd ih\^{a) is PT-symmetric, which sustains our previous assertion that the requirement 
of pseudo-Hermiticity covers a larger class of Hamiltonians, which have a positive spectrum. 
Once more we note that the dependence of h^^{a,g), H^j^{a,g) and ri{g) on the coupling 
constant g respects the identities ( 3.17| ). Notice that if the wave function cj) is vanishing at 



±00 this is no longer the case for ^ for odd m,. 

3.3.3 Harmonic oscillator perturbed by igx"^ 

In the examples discussed so far we were always able to construct explicitly the similarity 
transformation rj. However, when we start with a given non-Hermitian Hamiltonian this 
is not always possible. For instance, considering the simplest non-Hermitian perturbation 



of the harmonic oscillator by hi = gx |37, 44, 49] and their generalizations 



H^O^^^g^^^p2^^^2^-g^n (338) 

one has to resort at present to perturbation theory in order to construct rj. We adopt now 



the notation Sm,n from [^, 37, 44] for the totally symmetric polynomial in the m operators 
p and n operators x 

In the last expression we take the sum over the entire permutation group vr. Since the 
variables x and p are non-commutative this sum produces {m + n)\/ml/n\ non-equivalent 
terms. The first sum is the much simpler Weyl ordered version of this polynomial. Taking 
now in ( 3.21| ) the harmonic oscillator /i2(«) and h^ = gx'" as our starting point we solve 



Jn) 



( 3.21 ) for ql and compute to first order in perturbation theory 



J3) 2 /„ 2 



(9) 2 /^ 8 „ 48 ^ 64 „ 128 



'^i -al^^'^ + 3^^^'^+15^^^'^ + 35^^^'^ + 315^^^'°'' ^^'^^^ 



[("+i)/2] ^ r (A: - i - ^) 

<ii'^ = -V^ E (-a)^r(A. + i)r(i-§)^^'=-^-"+^-^^- ^^-^^^ 



10 
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We extrapolated here to the closed formula for all values of n, which we have verified up to 



n 



^ ' f~\ i-VT- /"vz-vci TTnt-l^ ■rl^/~\ c^ f~\\ ti r i f-\iri \-/~\-i iin rA ivi -(/ I [-i*~v t/^vtvi <-\it^it^i-v /i^ ' 



20. The expression for q\ agrees with the solution found in |37]. The remaining ql 



for n > 3 do not seem to be known in the literature. It is straightforward, but labourous 
to continue the analysis to higher orders. To next order we compute 



(3) 



^ (lE^^^'" + ^^^'^ + J^^^'^ - i^^'°) ' (^-^^^ 



The expression for q^ agrees precisely with the one found in |^, ^. Once again the 



expressions for higher values of n for gg seem to be unknown. From (3.15) we then 
compute the Hermitian counterpart to (p.38|) as 



h^io^.g) = Ip" + ?-^ - ^'^K.qt^] + 0{g'). (3.46) 



T 2 4 



The only commutator one needs to evaluate this is |52] 



\ — ^ i 1 77-1 T\ 

[x^,SrA = i 1. (_4)fc (2k + l)\{n-2k + l)\{r-2k + ^ySr-2k-i,s+n-2k~i. (3.47) 
where the upper limit of the sum is A(n, r) = min([(n + l)/2], [(r + l)/2]). We then obtain 

hn°ia,g) = -p^ + ^x^ + ^ Y^ Y,ct''S2ik-p-i),2in-P-k) (3-48) 

fc=l p=0 
with constants 

,.p = (_,.k^p+iVl^ r(2fc)r(n + i)r(fe-i-t) 

" ^ ^ 2 aHp r (fc + i) r (i - §) r {2k - 2p - 1) r (n - 2p) T{2p + 2) ^ ■ ' 

In particular, this reduces to 



^ 2^ +2" +2^ I 2'^^'^ + "'^°'^ "3 



h!^''i»,9) = oP' + o^' + o% 252,2 + a5o,4 - - (3.50) 



hf'^ia,g) = Ip^ + ^^^ + l^(t + "'-^cs - 4a5o,4 + 4a52,6 - 1652,2 + ^^4,4 ) (3.51) 



2 2 2 a'^ \ 5 3 

The expression for h^ {a,g) recovers the one found already in [p7| , ^ ]. This equivalence 
system has been studied extensively and here we shall elaborate further on it taking it as 
our prime example in the next section. Specifying now a = 1 the eigenvalue problem for 
the non-Hermitian counterpart of h^^{l,g), namely H-^^{l,g) was solved in ||5^ up to 
order g'^. The energy eigenvalues were found to be 

1 o^ 

en = n + - + ^(30n2 + 30n + 11) + 0(5^). (3.52) 

2 o 



The quite lengthy expression for corresponding wavefunctions $n(a^) niay be found in p3| , 
see formulae (3.2), (3.3) and (3.6) therein. In the next section we would like to use the 
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wavefunction for the Hermitian Hamiltonian h^ {l,g) instead, which we can simply com- 
pute from ( |2.5| ). With the help oi 
a differential operator in x-space 



pute from ( |2.5| ) . With the help of the explicit expression for q\ we may express the t] as 



r^ = l+^g 0^3 - xd^x^ - g' (x^ + 2x'd, - 38^ + ^d^ - ^xd^, - ^x'd^ + h^^ • 

(3.53) 
A somewhat lengthy but straightforward computation then yields 



•ng-xV2 



<l>^ix) = 7?$„(x) = ^__- [Hnix) - g^Pnix) + Oig^)] , (3.54) 



where the Hn{x) are the n Hermite polynomials, 

Pn{x) = ^ {2Hn-i{x) - (8n - 4)F„_2(x) + (2n + 3)F„+2(x) - ^Hn+i{x)\ (3.55) 

and Hn-p{x) = n{n — l)(n — 2) . . . (n — p + l)//„_p(x). The 4>n{x) are orthonormal 
wavefunction, which solve the Schrodinger equation up to order g^. We observe that 
despite the fact that the Hermitian Hamiltonian h^ is more complicated than its non- 
Hermitian counterpart H^^ , this is no longer true for their corresponding wavefunctions 
as 0„(x) takes on a much simpler form than $„(x). 

4. Time evolution for non-Hermitian Hamiltonians 

Next we consider genuinely time-dependent Hamiltonians. There have been some previous 
investigations in this direction |^, P, |l^, 11, 12, 13, 14|, which, however, do not make use 



of pseudo-Hermiticity. In addition, in many of these studies the precise meaning of the 



physical set up remains unclear. For instance, in [10|, no explanation is given about the 
meaning of making the mass time-dependent etc. Here we wish to address a more clear 
cut physical problem, namely one of the classical questions concerning the behaviour of a 
quantum mechanical system coupled to an external electromagnetic field. In particular, 
we have in mind an atom in a time-dependent linearly polarized electric field E{t) in the 
dipole approximation of finite duration r. In the length gauge, see section 4.2 for more 
discussions, this scenario is described by the Stark Hamiltonian and the time-dependent 
Schrodinger equation reads 



idtc^it) 



2 

y + y + xE{t) 



0(t) = [h + xE{t)] (f>{t) = hi{t)^it), (4.1) 



We follow here largely the notation of Q, ^, |5^. As the field is taken to be a pulse of 
finite duration we have h(j){0) = -E'(/>(0) and h(j){T) = E(J){t). With regard to our previous 
discussion we assume now that h has a non-Hermitian counterpart H which is in the same 
equivalence class, such that H = r]^^hrj. Hence this involves a potential for which we no 
longer demand that it is Hermitian, i.e. we allow V^ ^ V. Consequently also the resulting 
Stark Hamiltonian is non-Hermitian Hi{t) ^ Hj{t). 
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The central quantity of interest in this context is the time-evolution operator 

u{t,t') =Texp(-i dsh{s)] , (4.2) 

which evolves a wavefunction from a time t' to t, that is (j){t) = u{t,t')(f){t'). T denotes the 
time ordering. When h{s) is a self-adjoint operator in some Hilbert space, u{t, t') satisfies 
the relations 

idtu{t,t') = h(t)u{t,t'), u{t,t')u{t'X) = u{t,t") and u{t,t)=l. (4.3) 

Taking instead a Hamiltonian H{t) which is not self-adjoint and therefore its matrix repre- 
sentation is non-Hermitian these relations no longer hold. However, as we now demonstrate 
when H{t) is pseudo-Hermitian there is a simple modification of them. Acting adjointly 



with the time-independent operator rj on (4.3) and assuming that the similarity transfor- 



mation h = rjHrj ^ extends from the time-independent to the time-dependent system 

h{t) = vH{t)ri-\ (4.4) 

simply yields 

idtU{t, t') = H{t)U{t, t'), U{t, t')U{t', t") = U{t, t") and U{t, t) = I , (4.5) 

where we introduced the new time evolution operator U{t, t') associated to the non- 
Hermitian Hamiltonian H{t) as 

U(t,t') = T]-^uit,t')rj. (4.6) 

This time-evolution operator is quasi-pseudo-Hermitian 

UHt,t')=v'U-\t,t')ri-^ (4.7) 

which follows directly from w(t,t') = u~^{t,t'). The non-Hermitian counterpart H to the 
Hermitian Hamiltonian hi(t) as defined in (|4.1| ) results therefore to 



Hi{t) = H + 7j-^xE{t)r]. (4.8) 



The central assumption is here the validity of the similarity transformation (^^), which 
makes the formalism for the treatment of the non-Hermitian problem fairly straightforward. 
Of course we could also try to solve the problem for the situation when the electric field is 
coupled directly to the non-Hermitian Hamiltonian H, that means we take it to be of the 
form 

Hi{t)=H + xE{t). (4.9) 

In some special cases, namely when r]~^xE{t)r] = xE{t), this version is equivalent to ([4.^), 
but in general we require a new kind of formalism for this type of situation as we have now 
lost the equivalence relation ([4.4|). 
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4.1 Equivalent time-dependent pairs of Hamiltonians 

Let us illustrate the above for the concrete examples discussed in section 3. 

The simplest example is the generalized time-dependent Swanson Hamiltonian, which 
in the length gauge is simply of the form 

H^:Ua,g,t) = Hl^{a,g)+xE{t). (4.10) 



In this case the formulations ( [4.8D and ( |4.9| ) coincide as r/ ^xE{t)r] = xE{t). Its time- 
dependent Hermitian counterpart is therefore simply given by 

h^^^{a,g,t) = hl^^{a,g) + xE{t). (4.11) 

Thus for the Swanson Hamiltonian one obtains the same result whether one couples 
the electric field to h or H. 

For the perturbed anharmonic oscillators, this relation does no longer hold, since the 
similarity transformation ( 3.29[) does not commute with x, but instead induced a complex 



shift in X ^ x + ig/2. From (|4.8| ) the time-dependent versions of the anharmonic oscillators 
result to 

H^'''\a, g, t) = H^^'ia, g) + xE{t) + igE{t)/2, (4.12) 

which due to the last term is evidently different from the version ([4.S|). The time-dependent 
version of its Hermitian counterpart ( 3.32| ) is 



hi'^ia^g.t) = hi'^ia^g) + xE{t), (4.13) 

and one has now entirely different systems when coupling the electric field to h or H. 

The expressions become more complicated when we have non-trivial commutators 
between x and rj, as for the perturbed harmonic oscillator H^^{a,g), for which we only 
know the similarity transformation perturbatively. In that case the time dependent version 
becomes 

i/„^0(a,5,i) =i^r(«,5)+i^(t)E W^4"H^) (4.14) 

and we have to terminate the infinite sum according to the desired order of precision in 
powers of g. Using the commutator [x, Sm,n\ = imSm-i,n, which results as a special case of 
the commutator (p.47|), the first order in g is easily computed with the generic expression 



.(1) 



for qk', see ( |3^ ), to 

'('^+^)/^' 0F r(fc-i-f) 
t, (-«)^r(fc-i)r(i-t) 



H^''ia,g,t)=Hf^'^{a,g)+E{t) 



'^ 2^ (-rv^k r(h-l\r(l-n\ ^2k-2,n+l-2k 



(4.15) 



In particular we have 

H^''{a,g,t) = Hi'0^a,g)+xE{t) + i^E{t) (x' + -pA , (4.16) 

H^'^ia,g,t) = H^''{a,g)+xE{t) + i^E{t) (x^ + ^8^,2 + ^/) ■ (4.17) 
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The next order is more challenging as it involves commutators between different types of 
symmetric polynomials Sm,n and Sr^s for which an expression can be found, however, in 



m 



c- 1 _ • I I "v" 3i^2(-l - 2fc, -m, -s; n-2k,r- 2k; 1) ^ 

i^rn,n,^r,s\-^n-r- ^^ T {n - 2k)T {r - 2k)T {2k + 2) ^m+r-2k-l,n+s~2k-U 

(4.18) 
where 3F2 is a hypergeometric function. We will not present such calculations here. 

At this point one may wonder about the PT-symmetry of the non-Hermitian Hamil- 
tonians involved. For instance the term igE{t)/2 is only PT-symmetric ii E[—t) = —E(t), 
which means that it depends on the explicit form of the laser pulse. Taking for instance 
a typical pulse for a laser field with frequency oj, amplitude Eq and Gaussian envelop- 
ing function f{t), that is of the form E{t) = Eq sm(ujt) f (t) , would result in a PT-invariant 
Hamiltonian. However, the perfectly legitimate replacement sin(a;t) by cos(ci;t) would break 
the PT-invariance. Recall that in this context the electric field is treated classically. A 
discussion of PT-invariance for a full quantum electro dynamic setting may be found in 



1 56, ^]. For the physical application in mind, PT-invariance is, however, not a relevant 
issue here, since the pulse is always chosen such that /i#(0) = E^{0) and h^{T) = E^(t) 
and the eigenvalue problem is therefore only important in the time-independent case. We 



treat the full solution of (4T), the consequences on the non-Hermitian counterpart and 



dressed states [58| elsewhere [£9| 



4.2 Gauge transformations for non-Hermitian Hamiltonian systems 

For various applications it is extremely useful to transform the system to a different gauge. 
For instance, when having weak fields the length gauge is suitable as it usually involves the 
electric field just as an additional term, which is very useful for perturbation theory. The 
Kramers-Henneberger gauge is most useful when one wishes to exploit the periodicity of the 
field in Floquet analysis, especially for high frequencies. We now want to demonstrate how 
gauge transformations may be used for non-Hermitian Hamiltonian systems. Replacing 
for this purpose the wavefunction cp in the time-dependent Schrodinger equation related to 
some Hermitian Hamiltonian h hy (j) = a(t)~^(/)', with a{t) being some unitary operator, 
one obtains the well known identity, see e.g. ||5J, g, 55 1 



idt(b' = h'{t)4)' = [a{t)h{t)a{t)-^ + idta{t)a{t)-^] </>'. (4.19) 

Due to the relation (j) = r]^ it is straightforward to see that the gauge transformation for 
the non-Hermitian system results to 

idt<^' = H'{t)^' = [A{t)H{t)A{t)-'^ + idtA{t)A{t)-^] $', (4.20) 



where the similarity transformation (|4.4j ) extends to the gauge fields as well as to the gauge 
transformed time-dependent Hamiltonians 

a{t) = r]A{t)r]-^ and h{t) = r]H{t)r]^^. (4.21) 
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Note that the gauge transformations A{t) guarantee that physical observables remain in- 



variant, when computed using the generahzed inner product (2.6). 

In the context of laser-matter interaction, there are standard gauge transformations, 
from the length to the velocity gauge and from the velocity to the Kramers-Henneberger 
gauge 

az^4t)=e*W^ and a^^i^^W = e^'^^e-^^^P, (4.22) 

respectively, involving the classical momentum transfer 6(t), the classical displacement c{t) 
and the classical energy transfer d{t), from the laser field to the system in question. Such 
quantities are defined as 

b{t) = I dsE{s), c{t) = [ dsb{s) and d{t) = \ I dsh{sf. (4.23) 

Jo Jo 2 Jo 

In the Hermitian case, the Hamiltonians in the length, velocity and Kramers-Henneberger 
gauge are related as 

hi{p,x)-xE{t) =h^{p + b{t),x) = hKH{p,x + c{t)). (4.24) 

Physically, in the length gauge, the coupling with the field can be understood as a laser- 
induced dipole moment. In the velocity gauge, such a coupling appears as a shift p — > p—h{t) 
in the canonical momentum, corresponding to the well-known minimal coupling procedure. 
Finally, in the Kramers-Henneberger gauge, there is a displacement x — > x — c{t) in the 
coordinate x, which can be interpreted as time-dependent binding potential ||6^. For their 
pseudo-Hermitian counterparts H, one has in general 

Hi{p, x) - rj-^xE{t)r] / H^{p + b{t),x) / Hkh{p, x + c{t)). (4.25) 



Note that when in r/ = e'^'^ the operator q is linear in x or p, the equalities hold in ( 4.25| ). 



Otherwise, the similarity transformation will not induce a simple shift and will mix terms 



in X and p (for concrete examples, see section 4.2.1). 

We will compute the Hamiltonians discussed in section 3 in the velocity and the 
Kramers-Henneberger gauges, starting from their length-gauge counterparts (Sec. 4.1). 
Thereby, there exist two ways to proceed: either one applies the gauge transformations 
a{t) to the Hermitian Hamiltonians /i(t), and obtains its non-Hermitian counterpart em- 
ploying the similarity transformation 77, or one applies the transformations A{t) to the 
pseudo-Hermitian Hamiltonians H{t) directly. 

4.2.1 The generalized Swanson Hamiltonian 

For the generalized versions of the Swanson Hamiltonian, the similarity transformation 
7/ only depends on x, and therefore commutes with the transformation ai^v{t) from the 
length to the velocity gauge. This implies that ai^y{t) = Ai^y(t), so that the Swanson 
Hamiltonian Hn'm{c(,g,t) together with its Hermitian counterpart in the velocity gauge 
are easy to compute 

H^fm{a,g,t) = ^{p- b{t)f + |x" - i|(px™-i + x^-^p - 2b{t)x^-^), (4.26) 
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and 

hlU<^. g,t) = -{p- b{t)f + |x" + -5^x2— 2, (4.27) 

respectively. 

The computation of the time-dependent Swanson Hamiltonian in the Kramers-Henneberger 
gauge is slightly more involved, since the gauge transformation ay^KH{t) no longer com- 
mutes with the similarity transformation rj. Hence, 

A^^KHit) = exp[-^x'^]e^'^We-^"WP e^pi-x""] (4.28) 

771 m 

In this case, the time-dependent Swanson Hamiltonian and its Hermitian counterpart are 
computed to 

i/f;^^(a,g,t) = ^ + ^{x,{t)r + ^{x,{t)f"^'' -z^{px^-' +x"^-'p)-^x'^-', (4.29) 



and 



/if;^^(a,g,t) = ^ + |(xe(t))" + ^9'(xe(t))'"-^ (4.30) 



respectively, with Xc{t) = x — c{t). Note that in this case the equalities in ( |4.25| ) do not 
hold as there are terms in ( 4.29| ) occurring for which the displacement c{t) is absent. Such 



type of terms result from a shift p ^ p — igx"^ ^ in the momentum, caused by rj. This can 
be seen explicitly when ( [4.29| ) is re-written as 

<r («,ff,t) = ^P-'3f~'^" + |(x,(t))" + ^(Xe(t))^™"^ (4.31) 

4.2.2 Perturbed anharmonic oscillators 

For perturbed anharmonic oscillators, the similarity transformation -q depends on p. Hence, 
it no longer commutes with ai^v{t), so that 

Ai^^{t) = e-^^e^^W^e^P. (4.32) 

On the other hand, such a transformation commutes with a^^KH^t), from the velocity 
to the Kramers-Henneberger gauge. Therefore, Ay-.Kiii't) = cLv-^KHit)- For the non- 
Hermitian Hamiltonians, we then obtain 

r ("+l) i 

H^ (a,g,t)- + -X +y 2^ [-^ [2m -ir^'--''' 

rn=l ^ ^ ^ ' 

[ (" + !) ] 

F;^°'^^(a,,,t) = I- + \x.itT + ^ E (^) (2 J_ >2m-i(x.(t))"+--. 



In this case, the equality sign in (|4.25D hold in analogy to their hermitian counterparts 
( 4.24 ), which is expected, since g is a linear function of p. 

Once more we see from this that the non-Hermitian formulation exhibits no advantage 
over the Hermitian one. Even when in the time independent case the Hamiltonian K is 
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simpler than its Hermitian counterpart /i, this is spoiled by the introduction of the electric 
field. Thus in such a scenario, we have in /i a complicated potential term, but simple 
dependence on the electric field, whereas in H we have a simple potential but a complicated 
dependence on the electric field. Alternatively, one could add the field directly to H and 
thus keep both terms simple, but then the similarity transformation, which is already fixed 
by the time-independent part, will introduce non-Hermitian terms in h. Apart from this 
we have seen above that simplicity in the Hamiltonians does not imply simplicity in their 
eigenfunctions (see section 3.3.3). 

4.3 Perturbation theory for non-Hermitian Hamiltonian systems 

In most realistic situations, the time evolution of a physical system cannot be computed 
exactly. For instance, even in a Hermitian framework, it is in general not possible to solve 
the time-dependent Schrodinger equation for an atomic system with a binding potential 
V{x) subject to an external laser field E{t). Under these circumstances it is necessary to 
address the problem perturbatively. In this section, we will show how perturbation theory 
can be extended to a non-Hermitian framework. As a starting point, let us consider a 
time-dependent Hermitian Hamiltonian 

h{t) = ho{t) + hp{t) (4.33) 

where ho{t), hp{t) are also Hermitian and satisfy the time-dependent Schrodinger equation. 
Provided that the time-evolution operators associated to h{t) and /io(0 both satisfy the 
relation (14.31) , the time evolution operator u{t,t') associated to h can then be expressed by 
means of Du Hamel's formula [31, 54, g, ^^ 



Uh{t,t') = uo{t,t')-i Uh{t,s)hp{s)uo{s,t')ds. (4.34) 

Jt' 

By iterating Du Hamel's formula, one obtains a perturbative expansion for the time evolu- 
tion operator Uh{t, t') in hp <^ Hq. For instance, for a Hamiltonian of an atom in a potential 
V in the presence of an external laser field 

2 

h{t) = ^ + V{x)+xE{t), (4.35) 

one chooses hp{t) = xE{t) and hp = V{x) in the strong and weak field regime, respectively. 
In the latter case /io(i) is the Gordon- Volkov Hamiltonian, i.e., the Hamiltonian of a particle 



in the presence of the laser field only [62, S3[. 



As we have argued above, relations of the type ( |4.3| ) also hold for the time-evolution 
operator U{t,t') in (^^). Provided we can sensibly separate H{t) = H^it) + Hp{t) such 
that Uo{t,t') satisfies the relation ( [4. 5]) as well, Du Hamel's formula also holds for the 
non-Hermitian time-evolution operator 

Unit, t') = Uoit, t')-i [ UH{t, s)Hp{s)Uo{s, t')ds. (4.36) 

Jt' 
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The time-evolution operators may then be employed to compute various quantities of phys- 
ical interest, such as for instance the transition probability 



/ n 



($„|C/(t,0)|$, 



K</.Jn(t,0)|</.jr 



(4.37) 



from an eigenstate 



to \(j)^) of the Hermitian field- free Hamiltonian h or eigenstate 



|$m) to |$„) of the non-Hermitian field- free Hamiltonian H. We will consider first-order 



perturbation theory with respect to the external laser field amplitude £"0- Iterating ( 4.34| ) 
it follows that to this order the time-evolution operator can be approximated by 



'u(i)(f, 0) = no(t, 0) -i uo{t, s)xE{s)uo{s, 0)ds, 
Jo 

where uo(t, 0) = exp[— i/it]. The transition amplitude then reads 



{4'n\u{t,0)\(f)„ 



X p-is„t 
^ rt.rn '^ 



le 



-tent 



X 



Jo 



(4.38) 



(4.39) 



where the e„ are the eigenenergies. Form parity considerations it is evident that ((/>„| x \4>m) ¥" 
only when m, n are not both even or odd. 

Let us now employ these formulae for a concrete examples. 

4.4 Harmonic oscillator with a cubic non-Hermitian perturbation 

We consider the harmonic oscillator perturbed with a cubic non-Hermitian perturbation 
in the presence of a laser field 



1 



H^°'\l,g,t) =:H{t) = -p'^ + -x'^ + igx^ + r]-^xriE{t). 



(4.40) 



Up to order g^ this becomes with ( |3.40| ) 
Hit) 



p 



+ — +xE{t)+ig [x^ + x'^E{t) + 2p'^E{t)\+g'^E{t) [x^ - 2pxp\+0{g^) (4.41) 



As pointed out earlier, we observe that the additional term in ( 4.40| ), which contains the 
electric field has destroyed the simplicity of the time independent Hamiltonian. There is 
an additional problem with regard to perturbation theory, because we have lost the clear 
distinction of the potential term from the electric field term such that the separation into 
an -ffo(^) s-iid Hp{t) becomes more problematic, as now the two parameters g and Eq, which 
control the perturbative expansion, occur mixed, i.e. one has terms oc gEo, oc g'^Eo, etc. 

The Hermitian setting is much more clear cut and in addition the computations are 
far simpler as in the free field case the wavefunctions take on a simpler form as discussed 
in section 3.3.3. Thus when we consider the Hermitian counterpart of ( [4.41| ) instead 



h{t) =h + xE{t) = -p2 + -x2 + g^ 



|x* + 3S.,-i 



+ xE{t) 



(4.42) 



one is able to overcome this problem. We may now evaluate the transition probability 
( 4.37 ) for the eigenstates \(t>m)^ l'?^n.) of the field-free Hamiltonian hi^^{\,g) see ( 3.50| ), up 
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to second order in g, such that uo(i, 0) = exp[— z/ig t]. Choosing next a concrete form for 
the laser field E{t) = EQSiii{ijjt), that is a monochromatic driving field of frequency u) and 



field amplitude Eq, we may compute (4.39) for the Hamiltonian (4.42). Our results are 
presented in figure 1. 
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Figure 1: (Color online) Transition probability for a perturbed and unperturbed harmonic oscillator 
in the presence of a monochromatic laser field, as functions of the field frequency (panel (a)), and 
of the time (panels (b) and (c)). The perturbed and the unperturbed labels refer to expansions up 
to second and zero-th order in the parameter g respectively. We consider the transition from the 
energy level n = 15 to m = 16 to first-order perturbation theory with respect to the external laser 
field. The field amplitude is taken to be Eq = 0.003 a.u. and the coupling constant is chosen as 
g — 0.04. The pulse length r and the frequency cu are indicated in the figure. 



Panel (a) displays the transition probabilities |i;^i5) -^ l^Pw) or (|$i5) — > l^ie)) up 
to first order perturbation theory in Eq, when the system is subjected to a pulse of con- 
stant duration r, but with varying field frequency. For comparison, we also consider the 
unperturbed harmonic oscillator /iq in the presence of an external laser field, i.e. h{t) in 
0. Our choice of relatively highly excited states is motivated by the fact that 



(I4.42D fore/ 

according to ( |3.52 ) the difference between the perturbed and unperturbed system should 

be more pronounced for larger values of n. 

We expect to find that the system absorbs a single photon of frequency tu = Sn+i — £m 
in order to make a transition from the initial state |(/>„) to the final state If/'^+i). Thus for 
the unpertubed harmonic oscillator we expect a peak at u;(/io) = £«+! — £« = 1 and for the 



perturbed system we find from (3.52) that the peak should be at 



15 
Lo{h^^,n,g) = En+i -£„ = ! + 9^ y(™ + l)/2. 



(4.43) 



We evaluate from this io{h^ , 15,0.04) = 1.192, which agrees with our numerical calcula- 
tion of the expression ( [4.39D , resulting to 1.190. 
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Next we fix the frequency of the laser field to be resonant with the transition frequency, 
but vary the duration of the pulse. Prom standard computations (see e.g. [^]) it follows 
that the resonance probability should increase with t^ when the system is tuned to the 
transition frequency. This behaviour is confirmed by our perturbative calculations. In 
Panel (b) of figure 1 one clearly observes that that for the perturbed harmonic oscillator, 
the transition probability increases approximately quadratically in time, whereas the un- 
perturbed system does not exhibit this behaviour as it is off-resonance. In panel (c) the 
roles of the two systems are exchanged and we are now at the resonance frequency of the 
unperturbed system whereas the perturbed system is off-resonance. 



5. Conclusions 

We have constructed various new equivalence pairs, which relate non-Hermitian Hamil- 
tonians to their Hermitian counterparts. Our construction scheme is general and can be 
employed to compute further pairs not considered this far. 

We have demonstrated that when demanding the same similarity transformation to 
hold for the time-dependent Hamiltonian system and their time-independent counterpart, 
it is straightforward to develop a framework which describes the time evolution for non- 
Hermitian Hamiltonian systems. Despite the possibility to compute all relevant quantities 
in the non-Hermitian framework it turned out that it is usually easier to resort to the 
equivalent Hermitian formulation of the same systems and perform the evaluations in that 
context. In the future, it would be very interesting to investigate systems which constitute 
equivalence pairs in the time-independent case, but have the equivalence relation broken 
in the time-dependent scenario |5^, such as (|4.9|) . 

As was already remarked by various authors before, one may question the usefulness 
of PT-symmetry altogether. Pirst, despite the fact that it is a symmetry, it is does not 
guarantee that the spectrum of the non-Hermitian Hamiltonian will be real, due to the anti- 
linear nature of the PT-operator. Second, in the end it will come down to studying pseudo- 
Hermitian Hamiltonians, which not only constitute a wider class of systems, but in addition 
do not suffer from the shortcoming that the spectrum might not be positive and real after 
all. Third, we have used the fact that the time-independent non-Hermitian Hamiltonian is 
pseudo-Hermitian in formulating the time-evolution operator U . Just having PT-symmetry 
as the only principle at one's disposal would make this task very difficult. Of course using 
the relation rf = PC one may re-express all quantities in terms of CPT-operators, but 
that would really mean to use the similarity transformations as a construction principle. 
Pourth, apart from the time-independent Hamiltonian all quantities seem to be simpler 
in the Hermitian formulation, e.g. see above for the wavefunctions, the time-dependent 
Hamiltonians in their various gauges and perturbation theory. 

Despite its limited constructive power PT-symmetry remains useful in the sense that 
it is a very simple and transparent property, which can be read off directly from the 
Hamiltonian and thus constitutes a tool which can be used to identify potentially interesting 
non-Hermitian Hamiltonians. 
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